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Abstract. We present observations of four nearby early-type galaxies with previously known nuclear stellar disks 
using two instruments on-board the Hubble Space Telescope. We observed NGC 4128, NGC 4612, and NGC 5308 
with the Wide Field Planetary Camera 2, and the same three galaxies, plus NGC 4570, with the Space Telescope 
Imaging Spectrograph. We have detected a red nucleus in NGC 4128, a blue nucleus in NGC 4621, and a blue 
disk in NGC 5308. Additionally, we have discovered a blue disk-like feature with position angle ~ 15° from the 
major axis in NGC 4621. In NGC 5308 there is evidence for a blue region along the minor axis. We discovered a 
blue transient on the images of NGC 4128 at position 0'.'14 west and 0'.'32 north from the nucleus. The extracted 
kinematic profiles belong to two groups: fast (NGC 4570 and NGC 5308) and kinematically disturbed rotators 
(NGC 4128 and NGC 4621). We report the discovery of a kinematically decoupled core in NGC 4128. Galaxies 
have mostly old (10-14 Gyr) stellar populations with large spread in metallicities (sub- to super-solar). We discuss 
the possible formation scenarios, including bar-driven secular evolution and the influence of mergers, which can 
explain the observed color and kinematic features. 
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1. Introduction 

The decade of Hubble Space Telescope (HST) observa- 
tions have revealed the existence of small scale nuclear 
stellar disks in early-type galaxies. This discovery was 
an important step in the long process of recognizing the 
complexity of early-type galaxies. Ground-based studies 
^ , preceding the HST era, having lower resolution and pol- 
luted by typically > 1" seeing, alre ady recognized two 
, distinct classes of elliptical galaxies l|Davies et al.lll983|) 
that differed in photometric appearance - disky vs. boxy 
- and kinem atic propertie s - rotationally vs. pressure 
supported - l)Bendeil llQS'S iBender et aPligSgl) . Follow- 
up studies discovered the existence of embedded stellar 
disks in elliptical and lenticul ar galaxies ijScorza fc Benderl 
Il995t ISeifert fc Scorzalfl996(l . These disks, although sim- 
ilar to their counterparts in spiral and SO galaxies, have 
smaller scale length and higher central surface brightness. 
They often do not follow the same exponential profiles, 
and are closer to r^/^ profiles, refiecting formation in a 
different potential: dark halo potentials for disks in late- 
types, and bulge potenti als for disks in early-type galaxies 
llScorza fc Bended Il995h . The existence of the embedded 
disks also supports the idea of the morphological con- 



nection between spira l, lenticular and elliptical galaxies 
I Kormendv fc Bended ll996) . 



The properties of the nuclear regions (inner few 
100 pc) of early-type galaxies are, however, not eas- 
ily accessible from the ground. High resolution imag- 
ing surveys with HST discovered s mall scale nuclear 
stellar disks in ear ly type galaxies j jaffc et a l.l Il994t 

n 



van den Bosch et all Il994[ iLaiier et al.ni995,: .Rest et alJ 
200 ll hereafter ROll. They were followed by detailed 



photometric and kinematic stu dies on a few individual 
objec ts , principally NGC 4342 JScorza fc van den BoschI 
Il998t Ivan den Bosch et alJ 119981 hereafter B J M98)T 
NGC 4570 (BJM98; 'Scorz a fc van d en Bosch* 1998; 
van den Bosch fc EmseUem..l998ll. NGC 459 4 (Burkhead 
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198fit lKormendvlll98?^ lEmsellem et alJll99fil) NGC 7332 
l|Fisher et al.ll994LlFalc6n-Barroso et alJ2004|) . A detailed 
study of earl y-type galaxies with k inematicaly distinct 
components |Carollo et al.l Il997albl) found photometric 
evidences for faint nuclear stellar disks in a number of 
dust free galaxies. 

The next step was a search for embedded nuclear stel- 
lar disks in bulges of spiral galaxies. The high resolution 
studies of spiral galaxies with HST showed that a sig- 
nificant fraction of galaxies classified as early-type spi- 
rals have a rich variety of central properties, and show 
little evidence for r^ /*^ law expect ed for smooth bulges 
i Carollo et af] |l99i). Similarly, iBalcells et all ll2003^ 
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found moderately large fraction (34%) of nuclear bars or 
disks in their HST near-infrared survey of SO - She ga lax- 
ies. Using archival HST imaging. IPizzella et al.l l|2002fl re- 
ported evidence for nuclear disks in three early-type spi- 
rals and concluded that the disks are restricted to SO and 
unbarred spiral galaxies. Having in mind that nuclear stel- 
lar disks are detect able only when seen nearly edge on 
ijRix fc Whitell990(l . they appear to be very common, per- 
haps universal, in flattened ellipticals and SOs. 

The large fraction of detected nuclear disks in early- 
type galaxies (51% in the ROl sample) presents the ques- 
tions: how and when did the nuclear stellar disks form? 
Nuclear stellar disks are found in SO and disky ellipti- 
cals, but they are not simple extensions of large scale 
disks to the center of the galaxy. Often there is clear pho- 
tomet ric and kinematic evidences for double disk struc- 
tures llvan den Bosch et alJll994t IScorza fc Bended Il995l 
IScorza jT^iTderr^oschf Il998|) . where the double disk 
structures are represented by two morphologically sepa- 
rated disks, having different scale lengths, lying in nearly 
same plane, but possibly with different inclinations, and 
having an inner/ou ter separation radius be tween the disks. 
On the other hand. lSrwin fc SDarkelll2002l) show evidences 
that some inner disks seen in edge-on galaxies could be 
bars mistaken for disks. In any case, the two dynamically 
different structures are not easily distinguished in all cases. 

Inner disks are also found inside bars or rings 
llvan den Bosch fc Emsellemlll998t lErwin fc Sparkelll999t 
lErwin et al.ll2003(l . This is important for understanding 
their formation. Nuclear stellar disks could be the result 
of mergers in hierarchical galaxy formation scenario: ac- 
cretion of gas during the merger which settles in the prin- 
cipal plane of the galaxy and then makes stars. On the 
other hand, disks could be formed from the galaxy ma- 
terial transported to the nucleus by a bar, or perhaps 
from a mixture of these processes, in which a bar fuels 
the center effectively with gas captured at some previous 
epoch. Whatever scenario we choose, it has to be consis- 
tent with the high metallicities seen in the disks (BJM98; 
lEmsellem et al... IPOS') as well as their blue colors (BJM98; 
Kormendy et al. 2002) implying younger stellar popula- 
tions. 

Nuclear stellar disks with their cold dynamical 
properties and high surface brightness provide an 
excellent measure of the central mass-to-light ratio, 
as well as of the mass of the central black hole 



I van den Bosch fc de Zeeuwl I199(tI) . A few studies used 
this to determine the mass of the back holes in galax- 



ies with nuclear stellar disks 1 


<!ormendv et al.l 


.1996a.b 


ICretton fc van den Boschlll999l; 


Emsellem et al.lll999l). 



Early-type galaxies are also interesting for studying 
stellar populations. The absence of strong and continu- 
ous star-formation as well as emission line gas makes it 
easier to investigate the formation history and the con- 
nection between the photometric morphology, dynamical 
structures and corresponding stellar populations. 

In order to increase the available dataset and to in- 
vestigate the dependencies between the kinematics and 



Table 1. The properties of sample galaxies 



galaxy 


type 


Mb 


B-V 


Vrad 


PA 


D 


scale 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 4128 


5*0 


-19.89 


1.02 


2610 


58 


36.3 


175.7 


NGC 4570 


SO 


-20.39 


0.97 


1811 


159 


25.2 


121.9 


NGC 4621 


E 


-20.49 


0.97 


524 


165 


7.3 


35.3 


NGC 5308 


E-SO 


-20.38 


0.93 


2299 


60 


31.9 


154.8 



Notes - Col. (1): galaxy name; Col. (2): morphological type; 
Col. (3): absolute B-magnitude; Col. (4): apparent B-V color 
within the effective aperture in which half of the B-flux is emit- 
ted; Col. (5): radial velocity (cz) in km s~^ corrected for LG 
infall onto Virgo; Col. (6): major axis position angle in de- 
grees; Col. (7): distance in Mpc, as derived from radial veloc- 
ity (column 5) using H ubble constant Ho=72 km s"'^ Mpc~^ 
llFreedman et alJl200l[l : Col. (8): distance scale in pc arcsec"^. 
Values listed in columns 2-6 are taken from Lyon/Meudon 
Extragalactic Database (LEDA). 



line-strengths, as well as to determine the mass of black 
holes, we obtained high resolution spectra of four galaxies 
known to have nuclear stellar disks from the ROl sam- 
ple: NGC 4128, NGC 4570, NGC 4621 and NGC 5308. 
In addition, we also imaged three of the galaxies, ex- 
cept NGC 4570 which was th oroughly investigated in 
the previous studies (BJM98. IScorza fc van den BoschI 
Il998t lvan den Bosch fc Emsellemlll998l) . In this paper we 
present data observed with two instruments on-board HST 
during Cycle 9 (Program ID 8667) and concentrate on the 
photometric and spectroscopic properties and dependen- 
cies. Dynamical modeling of the galaxies with the purpose 
of determining the masses of the central black holes will 
follow in a separate paper. 

Section |21 presents broad band photometry, data re- 
duction, isophotal analysis and color images. Section |21 
deals with spectroscopic observations, data reduction, ex- 
traction of kinematics and measurements of line-strengths. 
Section 0] presents a discussion of the results for individ- 
ual galaxies. The conclusions and summary of the work 
are presented in Section |31 

2. WFPC2 broad band imaging 

The Wide Field Planetary Camera 2 (WFPC2) observa- 
tions (Biretta et al. 1996) included imaging in V (F555W) 
and B band (F450W) . The general properties of the sam- 
ple galaxies are presented in Table ^ the observations are 
summarized in Table[21and the details of the filter proper- 
ties are listed in Table |21 The centers of the galaxies were 
positioned on the PC CCD. The size of the PC CCD is 
800 X 800 pixels of 0'.'0455 x 0'.'0455. All exposures were 
taken with the telescope in fine lock. In addition to newly 
acquired data, we used existing archival I (F814W) band 
images for NGC 4621 (Program ID 8212, PI Ajhar) and 
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Table 2. The summary of HST/WFPC2 observations 



galaxy filter date time [s] # exp 



NGC 4128 F450W 14.05.2001 1400 2 

F555W 14.05.2001 800 2 

NGC 4621 F450W 14.05.2001 1200 2 

F555W 14.05.2001 800 2 

NGC 5308 F450W 07.05.2001 1400 2 

F555W 07.05.2001 800 2 



Table 3. The properties of the WFPC2 broad band fihers 
used (from WFPC2 Instrument Handbook), showing the 
central wavelength and the size of the wavelength window 
covered by the given filters 



filter Acen[A] AA[A] band 

F450W 4410 925 B 
F555W 5202 1223 V 



NGC 5308 (Program ID 5512, PI Faber). There were no I 
band observations for NGC 4128 in the archive. 

2.1. Data reduction 

The images were reduced through the standard 
HST/WFPC2 pipehne. Upon request of the data, the On- 
The-Fly reprocessing system re-reduced the data using 
the best calibration files. The standard reduction steps 
include correction for analog to digital conversion error, 
bias and dark curre nt subtraction and fla t-fielding (for 
more description see zman et all 1)19951) 1. Our obser- 
vations were divided (CR-SPLIT) in two (per filter). To 
combine them and to remove the cosmic rays we used 
a set of IDL routines from The IDL Astronomy User's 
Library fc andsmanl Il993|) ^ The WFPC2 images (PC, 
WF2, WF3, and WF4 for both CR-SPLIT sections) are 
cross-correlated to determine a possible shift between the 
exposures, aligned and combined removing the cosmic rays 
using an IDL equivalent of the IRAF task CRREJ. 

We then constructed color images: B-I, V-I and B-V. 
To construct the color images we had to align the individ- 
ual images very precisely. This was achieved by rotating 
the original images for the difference in the telescope ori- 
entation angle, sub-sampling pixels by a factor of six and 
cross-correlating images to find the shift. After all shifts 
were applied, the images were rebinned to the original 
pixel size. Both images used for the construction of a color 
image were initially convolved with the PSF of the other 
image. The PSFs were constructed using Tiny Tim soft- 
ware (Krist & Hook 2001). The raw counts of the images 
were converted into Johnson-Cousins B,V and I magni- 
tudes following the guidelines given by Holtzman et a l. 
(1995), using the zero points as given by Dolphin ('2000V 
and iterating the calibration until convergence. Note that 
the iteration was not performed for the B-I color images, 
because there were no published transformation for the 
F450W filter using B and I. We estimate our relative pho- 
tometric accuracy to be w 0.02 mag, while the absolute 
uncertainty is « 0.05 mag (and « 0.1 for B-I color). In 
Figs. Hand El we present WFPC2 observations, isophotal 



Wttp://idlastro. gsfc.nasa.gov 

\http:/ / www, noao. edu/ staff/ dolphin/wfpc2_ calib, 



analysis, broadband and color images of three observed 
galaxies. 

2.2. Isophotal analysis 

In order to investigate the disky structure of the galax- 
ies we used the IRAF task 'ellipse' to perform isophotal 
fits to the light distributions. We measured the ellipticity 
and position angle of the isophotes, as a f unction of ra- 
dius. T he method (for a full description see I.Tedrzeiewskil 
l|l987|) ') first fits elliptical isophotes to a Fourier expan- 
sion of first and second order terms. The next step in the 
method is to measure the higher order terms of the Fourier 
expansion. The pure ellipse is given by the first two order 
terms in the expansion. Any non zero values of the higher 
order ( > 2) terms means a de viation from the perfect 
ellipse. IPeletier et al.l CsH) and ICoudfrooii et al. 1 lll994^ 
found that cos39 terms (63) are sensitive to the presence 
of dust (as well as the difference between the higher or- 
der terms in different bands), while the cosAO terms (64) 
describe the shape by disting uishing boxy (b^ < 0) from 
disky (&4 > 0) galaxies fe.g. lLaueilll985t lBendeJll988|) . 
The isophotal paramete rs for three galaxies are sho wn in 
FigH Different studies l)van den Bosch et al.lll99l ROl) 
showed that although the isophote parameters can be fit- 
ted down to 0'.'03, they are not reliable and only values 
> 0'.'2 should be used for analysis. 

The galaxies were selected on the basis of having a nu- 
clear stellar disk and our isopohotal analysis agrees well 
with the ROl results. For all three galaxies the 63 terms are 
consistent with being zero in the reliable range (> 0'.'2). 
The 64 terms have clear disky deviations, but there are dif- 
ferences between the galaxies. NGC 4128 shows the small- 
est positive values of 64 coefficient, and is the only galaxy 
where 64 becomes negative. It drops below zero at a radius 
of ^ 3" from the center. The isophotes remain boxy in the 
rest of the investigated range (3 - 10"). The fourth order 
term in NGC 4621 is positive in the investigated range, al- 
though it starts to decrease beyond 3". NGC 5308 shows 
the sharpest rise of the 64 coefficient, but it drops to zero 
around 7" before it rises again. The analysis of photomet- 
ric higher order terms suggests the nuclear stellar disks are 
separated from the large scale disks. 

The coefficients of this isophotal analysis are one- 
dimensional representations of 2D structures in the galax- 
ies. In order to see more clearly the nuclear disks we 
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NGC 4128 NGC 4621 NGC 5308 
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-0.1 1 t I ^ i 

0.1 1.0 10.00.1 1.0 10.00.1 1.0 10.0 

r [arcsec] r [arcsec] r [arcsec] 

Fig. 1. Isophotal analysis results for NGC 4128, NGC 4621 and NGC 5308 in two observed filters: B (F450W) and V 
(F555W). First row: surface brightness profiles. Filled symbols correspond to B and open symbols to V filter. Errors 
are smaller then the symbols. Second row: position angles measured east of north; in the case of NGC 4621 we added 
90° to the measured position angle (east of north) for presentation purposes. Third row: ellipticity. Fourth row: higher 
order parameters 64 (the coefficient of cos(40)), showing deviations from perfect ellipses. 



constructed residual B and V images by subtracting a 
model galaxy from the original. The model galaxy was 
constructed using the IRAF task 'bmodel'. It had the 
same luminosity profile as the original galaxy, but it was 
constructed from perfect ellipses (first two terms in the 
Fourier expansion). The contour maps, shown in the mid- 
dle row in Fig. El reveal the disky structure that is respon- 



sible for the existence of the higher order terms. These 
disky structures do not necessarily show the full disks, 
because some information on the disks is also contained 
in the ellipticity, which was subtracted by the model. 
However, to the first approximation, these structures are 
a good representation of the relatively faint nuclear stel- 
lar disks in the observed galaxies. The disks have various 
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sizes between 2 - 5", corresponding to 150 ~ 600 pc, and 
are in a good agreement with the spatial extent of positive 
values of 64 terms. The noise structure in Fig. |21 perpen- 
dicular to the disks (along the minor axis) is not real and 
is an artifact of subtracting a perfect elliptical structure 
from a very disky one. 

2.3. Broad-band color images 

We constructed several broad band images using our data 
and archival images as described in Section 12.11 In the 
bottom two rows of Fig. |21we show B-V and B-I images 
for all three galaxies except for NGC 4128 which does not 
have archival I band images. These color images are color 
coded such that lighter shades indicate redder colors. 

2.3.1. B-V color images 

The nuclei of the galaxies show different structures in the 
B-V color. The nucleus of NGC 4128 is redder then the 
surrounding bulge, NGC 4621 has a convincingly blue nu- 
cleus, while NGC 5308 is rather uniform with a slightly 
bluer central pixel. The magnitudes measured on the con- 
structed color images are presented in Table 0] We mea- 
sured the color within a small circular aperture (6 pixels in 
diameter) to point out the subtle differences in colors be- 
tween different regions and features. The blue feature on 
the color image of NGC 4128 corresponds to a transient 
object detected in the galaxy. The nature of the transient 
is not known (see Appendix B for a detailed discussion). 
In NGC 4128 isochromes trace the nuclear disk, while in 
the case of NGC 5308 there is weak evidence of a blue 
disk outside O'.'S of the nucleus. The difference is « 0.02 
magnitudes. NGC 4621 is a special case different from 
both previous galaxies. Here we have a prominent blue 
nucleus on top of the red bulge. The other particularity of 
this galaxy is in the blue component that stretches south- 
wards from the nucleus. The average difference between 
this component and the rest of the bulge is again sa 0.02. 

2.3.2. B,V-I color images 

A bigger color difference is seen in the B-I and V-I images 
which we constructed for NGC 4621 and NGC 5308. The 
archival I-band images of NGC 4621 were of good quality, 
while images of NGC 5308 were saturated, and the very 
central parts of the images are not reliable. 

The interesting blue feature southwards from the nu- 
cleus on the B-I image of NGC 4621 is more prominent and 
we can quantify its extent. The color difference between 
the red bulge-like background and the blue patch is here 
« 0.12 mag in its extreme. The extent of the feature is 
« 1'.'35 on each side of the nucleus, being somewhat bluer 
on the south side. The feature lies generally in the north- 
south direction making an angle of 15° to the major axis. 
It is in no obvious plane of symmetry of the galaxy. The 
rest of the bulge is red, but a region of enhanced red color 



Table 4. Color in magnitudes measured within apertures 
of 6 pixels (0'.'273) in diameter 



galaxy 


B 


- V 


V - 


I 


B - 


I 




center 


average 


center 


BF 


bulge 


BF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 4128 


1.00 


0.95 










NGC 4621 


0.87 


0.99 


1.22 


1.16 


2.14 


2.05 


NGC 5308 


0.97 


0.94 











Notes - Col. (1): galaxy name; Col. (2): B-V color measured 
at the center of the galaxies; Col. (3): B-V color averaged over 
8 apertures placed around center on a square grid centered 
on the nucleus with size 2 x 14 pixels (~ l'.'27), except for 
NGC 4128 where the blue feature next to the center was ex- 
cluded; Col. (4): V-I color measured at the center of galaxy; 
Col. (5): V-I color measured at the blue feature 0'.'90 from the 
center; Col. (6): B-I color measured at bulge of the galaxy; 
Col. (7): B-I color measured at the blue feature 0'.'90 from the 
center. 



lies in the east-west direction. The angle between this red- 
dest region and the major axis of about 103°, making the 
angle between the bluest and the reddest regions about 
118°. T he V-I color image c onfirms this finding, also no- 
ticed bv lWernh et all l)2002|) . hereafter WEC02, on their 
V-I images. We postpone the discussion and interpretation 
of these features to Section 

NGC 5308 is an equally interesting case. The center 
of the I band image was saturated and the blue nucleus 
is an artifact. Generally, one has to be very careful in 
interpreting the saturated images. Aware of problems in 
dealing with saturated images, we used them to verify the 
faint suggestion from the B-V images of the blue disk. We 
do not consider the central 0'.'5 of the B-I and V-I color 
images (the area somewhat larger than the blue dip in the 
images), but concentrate on the larger-scale features. 

Along the minor axis of the galaxy, the B-I color im- 
age also reveals blue regions on each side of the nucleus, as 
well as a red feature in the east-west direction. These fea- 
tures were not anticipated and we checked whether they 
are real or artifacts of saturation since they are close to 
the nucleus (r < 1"). We first constructed color image 
from unconvolved F4550W and F814W images. On the 
resulting color image, there was a hint of the east-west 
red feature. The other approach included deconvolution 
of the B and V images using Richardson-Lucy algorithm 
with 20 iterations. The resulting images were convolved 
with the PSF of the I image, and used to make color im- 
ages (B-I and V-I). On both color images, next to the red 
east-west feature, we also detected the blue region on the 
minor axis. These tests suggest the color features on the 
last panel of the Fig. |3 are real, although could be aug- 
mented by the convolution process due to its proximity to 
the saturated nucleus. 
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NGC4128 NGC4621 NGC5308 




arcsec arcsec 
Fig. 2. WFPC2 observation of NGC 4128, NGC 4621, NGC 5308. First row: images in the V (F555W) filter. The 
arrows and their associated dashes mark the North and East orientation of the images. Second row: residual stellar disks 
after subtraction of perfect elliptical model galaxy (in V band). Contours were slightly smoothed and logarithmically 
scaled with steps of 0.5 magnitudes. Third row: B-V images, with grayscale such that brighter means redder. The blue 
feature next to the center of NGC 4128 is due to the transient and is discussed in Appendix B. The plotted range is 
(in magnitudes) 0.6 (black) to 1.0 (white) for NGC 4128, 0.5 (black) to 1.05 (white) for NGC 4621, and 0.5 (black) to 
0.98 (white) NGC 5308. Fourth row: B-I images. The central blue feature on the NGC 5308 color map is an artifact 
of saturation and was not considered in the analysis. Range is from 1.93 (black) to 2.16 (white) for NGC 4621, and 
1.1 (black) to 2.07 (white) for NGC 5308. 



The remaining and real (clearly visible on all tost 
images), thin blue feature follows the major axis of the 
galaxy, along which there is clear evidence for a very thin 
disk on all scales, from the nucleus outwards. The posi- 



tion angle of the thin blue component is the same as of 
the nuclear stellar disk. 

The marginal difference in color (« 0.01 — 0.02 mag) 
is almost razor sharp and it looks like the signature of the 
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Table 5. Summary of HST/STIS observations 






□rcsec 



Fig. 3. Comparison of B-I color profiles extracted (and 
smoothed) along the disk (thin line) and parallel to 
disk (thick line) in NGC 5308. The slit-like apertures 
were 200x1 pixels in size. One slit (thin line) was placed 
along the major axis (and the disk), while two other slits 
(averaged together and presented by the thick line) were 
placed parallel to the major axis, 10 pixels (0'.'45) above 
and below it. Vertical dashed lines show the nuclear region 
excluded from the measurement due to the saturation 
effects. 

thin disk visible also in the residual image (second row of 
Fig. 12). We compared the sizes of the disk in the residual 
image and the blue feature in the B-I color image. The 
comparison is made by extracting and averaging together 
several profiles of intensity and color perpendicular to the 
disk, on both sides of the nucleus (avoiding the central 
2"). The final profiles were fitted with Gaussians. The size 
(FWHM) of the disk feature on the residual disk image is 
~ 5.6 pixels and the size of the color feature is ~ 4.4 pixels. 
These numbers correspond to 0'.'25 and 0'.'20 respectively, 
and they are in good agreement, enforcing the connection 
between the components. At the distance of NGC 5308, 
the blue component in the color image is approximately 
30 pc thick. 

A way to quantify the relative difference in color be- 
tween the disk and the bulge is shown on Fig. |3| We 
measured the color along slit-like apertures (1 pixel wide) 
along the disk and parallel to it, on both sides of the nu- 
cleus. The central 1" of all slits were omitted, and the 
two slits, positioned on each side of the central slit, were 
averaged and presented as one color profile. The disk is 
clearly bluer then the bulge in the inner 8". Beyond 3'.'5 
on each side of the center colors of the disk and bulge 
become similar. 

Briefly summarizing, we list below the important ob- 
served color features to which we refer later in the text: 



galaxy 


Silt 




date 


time 


# exp 


A 

^cen 




(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


{<) 


NGC 4128 


ccn 


02 


12.2001 


2568 5 


3-1-3 


0.00 


— 112.9 




pos 


02 


12.2001 


2700 


3 


0.16 


-112.9 






02 


12.2001 


2697 


3 


—0.40 


-112.9 


NGC 4570 


cen 


01 


04.2001 


2249.5 


3-1-3 


-0.20 


152.1 




pos 


01 


04.2001 


2369 


2 


0.68 


152.1 




neg 


01 


04.2001 


2520 


2 


-0.56 


152.1 


NGC 4621 


cen 


01 


04.2001 


2260.5 


3-F3 


-0.12 


161.1 




pos 


01 


04.2001 


2380 


3 


0.40 


161.1 




neg 


01 


04.2001 


2520 


3 


-0.44 


161.1 


NGC 5308 


cen 


12 


07.2000 


2525.0 


3-f 3 


0.00 


60.4 




pos 


12 


07.2000 


2667 


3 


0.36 


60.4 




neg 


12 


07.2000 


2670 


3 


-0.36 


60.4 



Notes - Col. (1): galaxy name; Col. (2): the position of slit, 
cen - center, pos - positive offset, neg - negative offset (see 
text for details); Col. (3): date of observations; Col. (4): total 
exposure time in seconds (exposure time of central slits is the 
average time of all added observations); Col. (5); number of 
used observations; Col. (6): distance of the center of the slit 
from the galaxy nucleus in arcsec; Col (7): position angle of 
the slit in degrees east from north. 



NGC 4128 has (i) a red nucleus and (ii) a blue feature 
0'.'14 west and 0'.'32 north of the galaxy center; 

NGC 4621 has (i) a blue nucleus, (ii) an extended blue 
component with PA « 150° (east of north), and (iii) an 
extended red component with PA « 268° (east of north); 

NGC 5308 has (i) an extended blue component along the 
major axis and (ii) a wide blue component along the minor 
axis, (iii) an extended red component with an east-west 
orientation. 



3. STIS spectroscopy 

Spectra of all four galaxies were obtained using the Space 
Telescope Imaging Spectrograph (STIS) \ yith aperture 
52"xy 2 and grating G430M (for details see lKimble et alJ 
l|l998|) '). The observations are summarized in Table The 
configuration of STIS and the properties of the grating are 
listed in Table El With this setup we chose to observe Mg 
lines at 5180 A rather then Ca lines at 8700 A which 
are commonly used for extracting stellar kinematics. The 
reasons were: (i) the Mg lines provide simultaneous kine- 
matic data and a commonly used index of stellar mctal- 
licity, (ii) STIS shows problems with scattered light inside 
the CCD chip in the near infra-red which generates ar- 
tificial "wings" on spectral features, and (iii) the spatial 
resolution along the slit in V-band is better by almost a 
factor of two due to the decreased Airy disk size. 
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Table 6. The configuration of STIS and the properties of 
the grating 



Quantity 



Value 



Aperture 
Grating 
A-range [A] 

Acen [A] 

Scale AA [A pixeP^] 
Spatial scale [arcsec pixel" ^] 
Comparison line FWHM [pixel] 
R =A/AA 

Instrumental dispersion [km s^^] 



52x0.2 
G430M 
5050.4 - 5381.6 
5216 
0.28 
0.05 
2.9 
6461 
19.76 



3.1. Data reduction 

The galaxies were observed in a similar manner with a 
total of four orbits per galaxy. Spectra were taken at three 
parallel positions per galaxy during the four orbits. Each 
orbit was divided (CR-SPLIT) into 3 shorter exposures. 
Two orbits were used for the slit placed on the center of the 
galaxies (cen), along the major axis. Between orbits, the 
galaxy was shifted along the slit for about 0'.'2, or 4 pixels, 
to get a better estimate of detector sensitivity variations 
and to identify hot pixels. This strategy was not successful 
for the case of NGC 4128, where the measured shift was 
~ 1 pixel. The remaining two orbits were split between 
two slit positions on either side of the central slit, covering 
the bulge parallel to the nuclear stellar disk. One slit was 
targeted at the position +0'.'3 away from the central slit 
(positive offset - pos), and the other at the position -O'.'B 
away from the central slit (negative offset - neg). 

Most of the data reduction was performed by the 
HST/STIS calibration pipeline CALSTIS (Hodge et al. 
1998), including subtraction of overscan, bias and dark, 
then flat-fielding, hot pixel and cosmic-ray removal, ab- 
solute sensitivity calibration and wavelength calibration. 
CR-SPLIT data sets were combined automatically in the 
pipeline. The additional combination of the two central slit 
exposures was performed manually outside the pipeline in 
IRAF, using the STSDAS task MSCOMBINE. This task 
averages the exposures scaled by their exposure times, and 
combines the separate exposures using a robust sigma clip- 
ping rejection method. This was done on the files that 
have _crj extension, i.e., after pipeline co-addition of CR- 
SPLIT images and before the calibrations. The combined 
files were then returned to CALSTIS for the calibrations. 
The same procedure, however, was not possible on the side 
slits which were taken during only one orbit. As there were 
still some cosmic rays left after the pipeline reduction, we 
used Lapl acian Cosmic Ray Id entification (LAcosmic) de- 
veloped bv lvan Dokk7m]|)200ll) to remove them. LAcosmic 
was also applied on the _crj files. The detection limit for 
the outliers was 3.5cr. To improve the quality of the images 
and remove additional negative pixels we tried a few tech- 
niques. For spectra with very low signal-to-noise it was 



0.0 
1 .0 







Fig. 4. Galaxy light profile along the STIS sht. The full 
spectral range was used. Thin dark line is the profile along 
the cen slit, thick gray lines are profiles along the side slits: 
dashed line along the pos slit and full line along the neg 
slit. The profiles were normalized to the maximum of the 
cen slit in order to emphasize the difference of intensities. 

possible to compare different exposures (e.g. the two side 
slits) and to recognize the same negative pixels and create 
a mask of them. For spectra with higher signal to noise 
this was not effective and a different approach was used. 
Using a boxcar filter we smoothed each LAcosmic-filtered 
image. These images were subtracted from the correspond- 
ing original LAcosmic-filtered images to emphasizes the 
outlying pixels. They were flagged creating a mask image. 
Masked pixels were interpolated using IRAF task FIXPIX. 
The resulting images were returned to CALSTIS and pro- 
cessed to the end of the pipeline. 

The final STIS light profiles are shown in Fig.0] A no- 
ticeable feature is the difference in the intensity of the side 
slits. If the centering and shifting process worked properly, 
as the light profiles of the galaxies are quite symmetric, it 
is expected that the side slits should have very similar pro- 
files. This is true for the case of NGC 5308 and NGC 4621, 
suggesting those slits were on similar but opposite posi- 
tions. However, the other two galaxies show significant 
deviations. It is therefore necessary to find the exact po- 
sitions of all slits. 

We checked the actual positions of the slits by com- 
paring the light profiles from the STIS spectra with the 
WFPC2 images. The width of the slit (dispersion direc- 
tion) is 0'.'2, which corresponds to 4 pixels on the STIS 
CCD. We sub-sampled the image such that the slit width 
projects to 5 pixels, in order to center the slits more 
correctly. Summing up along the dispersion axis (x-axis 
on the CCD) we created a STIS light profile. This was 
repeated for all slits. The F555W images, which were 
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Fig. 5. Plots of versus the positions of the shts for each 
galaxy. The line connecting the diamonds is the relative 
of the central slit. The line connecting asterisks is the 
total relative obtained by including the side slits in 
calculations. 

used in the comparison, were accordingly re-sampled. In 
each case we scanned the WFPC2 image by a combi- 
nation of the three slits, independently varying the dis- 
tances between the slits. The comparison of the STIS and 
WFPC2 profile was expressed by the relative {pro- 
file(WFPC2)/profile(STIS) - 1). In this process we as- 
sumed that the position angle of the telescope did not 
change between different slit positions. The resulting 
estimates are shown in Fig. [S] If the central slit is not 
on the nucleus, the is expected to have double min- 
ima and it is hard to distinguish which one is correct 
(NGC 4621 and NGC 4570 are clear examples). However, 
using the additional light profiles of the side slits tight- 
ens the constraints producing one clear minimum, which 
corresponds to the position of the central slit. The uncer- 
tainty of our estimate is 0'.'04. The positions of the central 
and the side slits with respect to the nucleus are given in 
last two columns of Table El The slits in NGC 5308 and 
NGC 4128 are centered on the galaxy nucleus, while for 
NGC 4621 and NGC 4570 the slits were offset. The side 
slits were roughly on the requested positions for NGC 4621 
and NGC 5308. In the case of NGC 4128, the side slits are 
the least symmetrically positioned and the pos slit is al- 
most coincident with the central slit; while in the case 
of the NGC 4570 the positions of the side slits are the 
farthest apart, as suggested from the profiles (Fig.0J). 

The general characteristic of the spectroscopic data 
is their low signal-to- noise ratio [S/N). Only the major 
axis spectra have an S/N sufficient for the extraction of 
kinematics and line-strengths as a function of radius. The 
side slits are much noisier and no kinematic measurements 



were possible. It was possible however to extract line- 
strength information from a few central rows, summed 
together to increase S/N creating one spectrum per side 
slit per galaxy. 

3.2. Stellar kinematics 

All available information about the stellar kinematic prop- 
erties of galaxies are given by the line-of-sight velocity 
distribution (LOSVD). The process of extracting kinemat- 
ics is based on the deconvolution of the observed galaxy 
spectra in order to recover the full LOSVD. The idea 
behind this is that the galaxy spectrum can be repro- 
duced using a combination of several representative stellar 
spectra convolved with the true LOSVD. Unfortunately, 
the LOSVD is not a priori known and the deconvolu- 
tion process is ill-determined, being heavily dependent 
on the quality of the data. Over the last thirty years a 
number of methods were invented to tackle the problem 
and deliver the b est possible estimates of the LOSVD (see 
Ide Bruvne et alJ (|2p03) for an overview of methods). Here 
we choose to use a parametric method operating in the 
pixel space because of the low S/N of our data and very 
short wavelength ran ge. We use the penalized pixe l fit- 
ting (pPXF) method ICappellari fc Ems ellcm"2004- We 
deriv e the LOSVD parameterized by a Gauss-Herm ite se- 
ries (jvan der Marel k Franxl Il993 ICerhardl Il993(l . The 
method finds the best fit to a galaxy spectrum by convolv- 
ing an optimal template spectrum with the corresponding 
LOSVD given by the mean velocity V and velocity disper- 
sion cr, as well as higher order Gauss-Hermite moments 
and /14. The higher order moments measure asymmetric 
and symmetric deviation of the LOSVD from a Gaussian 
respectively. 

An element which can heavily influence the extracted 
kinematics is the stellar template used to convolve the 
LOSVD to reproduce the galaxy spectra. T here are met h- 
ods, such as Fourier correlation quotient dBendeiUlQQfil) 
or Cross-correlation method (Statler 199^, which are less 
sensitive to template mismatch. Pixel fitting techniques 
are much more sensitive to template mismatch and it is 
crucial to have a good stellar template before starting the 
extraction. The usual way is to observe a number of repre- 
sentative stars (matching the spread in metallicity and age 
of stars in the observed galaxy) with the same instrumen- 
tal set-up and to build an optimal template as a weighted 
linear combination of the observed stellar spectra. 

After searching through the HST archive we decided 
to use a set of stellar population models instead of the 
one star from the archive that matched our set-up (the 
same grism being the most important, while size of the 
slit can be accounted for). Usin g singl e-metallicity stellar 
population models of IVazdekil l)l999l) we constructed a 
large stellar library from which to build the optimal stellar 
template. Each galaxy long-slit spectrum was summed up 
along the slit to make a higher S/N spectrum, which was 
used to obtain the optimal template. We also used additive 
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Fig. 6. Example of spectra for the four studied galaxies. 
From bottom up: NGC 4128, NGC 4570, NGC 4621 and 
NGC 5308. The spectra are shifted vertically to avoid over- 
lap. The thin black lines indicate spectra of galaxies. The 
dashed thick lines are broadened optimal templates. The 
dots below the spectra are residuals of the fit (difference 
between galaxy spectra and optimal templates). Vertical 
and horizontal lines shows the region used in the fit. Two 
solid vertical lines crossing the spectra show the spectral 
regions excluded from the fit which tested the influence of 
Mgb region to the extracted kinematics. 

Legendre polynomials to adapt the continuum shape of the 
templates. This optimal template was then used in the fit 
of the individual spectra along the slit. 

A disadvantage of the Vazdekis models is that they 
are of lower resolution than the STIS data. The FWHM 
of the Vazdekis library is 1.8 A compared to 0.8 A from 
STIS. This requires a degradation of our data by ^1.6 
A. Although certain information is in this way lost and 
the STIS spectral resolution is degraded, the smoothing 
of the data helps in removing noise and the library en- 
sures the extracted kinematics do not suffer from an im- 
portant systematic template mismatch effect. Examples 



of constructed optimal templates are shown in Fig.|Sl The 
presented spectra are the sum of the central 10 rows, hav- 
ing a high S/N which is needed for properly estimating the 
optimal template. The overplotted dashed lines are our re- 
sulting optimal templates for the galaxies, convolved with 
the determined LOSVD of the galaxy. Typically a few (2- 
3) old-type stars from the Vazdekis library were selected 
by the fitting routine for the optimal template. The resid- 
uals between the galaxy and the optimal template spectra 
are shown below each spectrum. The presented optimal 
templates were used to extract the kinematics from the 
spatially binned spectra. An alternative way would be to 
construct the optimal template for each spatial bin and 
then use this to extract kinematics in the same bin. This 
method is important for galaxies with stellar populations 
changing between bins (~ 0'.'05), but in this case the low 
S/N of the individual spectra do not justify this approach. 

Table summarizes the details about the spatial bins 
used for the extraction of kinematics. They were chosen 
after some experimenting as a compromise between the 
S/N and the spatial resolution. The galaxies have differ- 
ent surface brightnesses and, since the exposure times were 
similar, a unique scheme for all galaxies was not useful. For 
each galaxy we assumed a target S/N and we binned ac- 
cordingly. Generally, the spectra become too noisy to mea- 
sure the kinematics beyond 1". In some cases, the central 
few rows of spectra have S/N high enough for extraction 
of the higher order terms of LOSVD, but in general the 
S/N is too low. Hence, we decide to confine the extrac- 
tion to only the first two moments (assuming a Gaussian 
shape for LOSVD): mean stellar velocity {V) and velocity 
dispersion (cr). 

Another element which can heavily influence the re- 
sults of the extra ction is speci f ic to t he spectral region of 
the observations. iBarth et alJ l|2002f l compared the kine- 
matics of a number of galaxies extracted in two spec- 
tral regions: one around Mg5 lines and the other around 
Ca triplet. They found that if the metallicities of the 
galaxies and template stars are not well matched then 
direct template-fitting results are improved if the Mg6 
lines themselves are excluded from the fit and the veloc- 
ity dispersion is determined from the surrounding weaker 
lines. For galaxies with high velocity dispersion this will 
be more important because of the correlati on between the 
veloc i ty dispersion and the [Mg/Fe] ratio ( Worth ev et alJ 
ll 9fl2t iTra.ger et a,1.L ,1C)98: Kuntschner et a,l. 2001.1 . which 
increases the strength of the Mg6 lines relativ e to the sur- 
rounding Fe lines. Following the suggestion of^Barth et alJ 
(|2002} we also extracted kinematics excluding from the fit 
the Mg6 lines (the excluded regions are shown on Fig. [Bias 
vertical lines crossing the spectra). When there are signif- 
icant differences between the two sets of extracted kine- 
matics we used the set obtained by excluding the Mg6 lines 
from the fit for the further analysis and interpretations. 

The errors were estimated using Monte-Carlo simula- 
tions. The LOSVD parameters were derived from 100 re- 
alizations of the input spectrum, where the value at each 
pixel is taken from a Gaussian distribution with the mean 
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Table 7. Bins for kinematic extraction 



galaxy 


bin r 


[arcsec] 


width [pix] 


ran 


ge 


S/N 


NGC 4128 


center 


0.00 


1 


599 - 


599 


15 




rl 


0.05 


1 


600 - 


600 


14 




r2 


0.15 


2 


601 - 


603 


15 




r3 


0.25 


3 


604- 


607 


13 




r4 


0.50 


6 


608 - 


612 


12 




r5 


1.25 


24 


613 - 


637 


12 




11 


—0.05 


1 


598 — 


598 


14 




12 


-0.10 


1 


596 - 


597 


12 




13 


-0.15 


2 


593 - 


595 


14 




14 


-0.30 


3 


589 — 


592 


13 




15 


-0.50 


5 


583 - 


588 


12 




16 


-0.95 


13 


569 - 


582 


12 


NGC 4570 


center 


0.00 


1 


599 - 


599 


22 




rl 


0.05 


1 


600 - 


600 


21 




r2 


0.10 


1 


601 - 


601 


18 




r3 


0.15 


2 


602 - 


604 


20 




r4 


0.30 


3 


605 - 


608 


20 




r5 


0.50 


5 


609 - 


614 


19 




r6 


0.80 


8 


615 — 


623 


19 




11 


—0.05 


1 


598 — 


598 


19 




12 


-0.15 


2 


595 - 


597 


21 




13 


-0.25 


3 


591 - 


594 


21 




14 


—0.45 


4 


586 — 


590 


19 




15 


-0.70 


6 


579 - 


585 


18 


NGC 4621 


center 


0.00 


1 


599 - 


599 


40 




rl 


0.05 


1 


600 - 


600 


34 




r2 


0.10 


1 


601 - 


601 


25 




r3 


0.20 


2 


602 - 


604 


26 




r4 


0.35 


4 


605 - 


609 


27 




r5 


0.65 


8 


610 - 


618 


26 




11 


-0.05 


1 


598 - 


598 


36 




12 


—0.10 


1 


597 — 


597 


28 




13 


—0.15 


2 


594 — 


596 


29 




14 


-0.30 


3 


590 - 


593 


26 




15 


-0.50 


6 


583 - 


589 


27 




16 


—0.95 


11 


571 — 


582 


25 


NGC 5308 


center 


0.00 


1 


599 - 


599 


27 




rl 


0.05 


1 


600 - 


600 


23 




r2 


0.10 


1 


601 - 


601 


17 




r3 


0.20 


2 


602 - 


604 


17 




r4 


0.30 


3 


605 - 


608 


15 




r5 


0.55 


6 


609 - 


615 


15 




11 


-0.05 


1 


598 - 


598 


23 




12 


-0.10 


1 


597- 


597 


17 




13 


-0.15 


2 


594- 


596 


18 




14 


-0.30 


3 


590 - 


593 


17 




15 


-0.50 


5 


584- 


589 


16 




16 


-0.90 


11 


572 - 


583 


15 



of the initial spectrum and standard deviation given by a 
robust sigma estimate of the residual of the fit to the initial 
spectrum. Fig.|B|shows an example of the residuals used to 
estimate the standard deviation used in Monte-Carlo sim- 
ulations (dots under the spectra) . All realizations provide 



a distribution of values from which we estimate the Icr 
confidence limits. The values of the extracted kinematics 
are presented in Tables A. 1-2 of Appendix A and shown 
in Figs.|Hl-ini 

All galaxies except NGC 4621 show rather fast major 
axis rotation. NGC 4621 is a special case with a previously 
discovered counter-rotation in the center (WEC02). There 
are some differences between the kinematics extracted fit- 
ting the full spectral range and excluding the Mg6 lines. 
They are the strongest for NGC 4621 and NGC 4570. The 
somewhat larger error bars of the kinematic measurements 
obtained not fitting the Mg6 region are the consequence 
of lowering the S/N by excluding the dominant spectral 
feature. We postpone detailed description of all kinematic 
curves to Sectional 

3.3. Line strengths 

The spectral range of our observations is very limited cov- 
ering only the Mg6 and Fe5270 Lick/IDS i ndices (for def- 
inition of Lick/IDS sys t em and indices se e'Burst ein et alJ 
(|1984^ : IWort.hev et a,lJ l)l994j) : iTrager et a,l.. The 
red continuum pass band of the Fe5270 index is trun- 
cated by the edge of our spectral range and this in- 
dex cannot be used in its defined form. A similar case 

is f ound in Kuntschner et al. 2004 (in preparation) 

and lFalcon-Barroso et al.l (2004) where Fe5270 cannot be 
mapped over the who le field-of-view of the integral-field 
spectrograph SAURON ( Bacon et al.ll200l]) due to the vary- 
ing bandpass of the SAURON instrument. In their case, 
Kuntschner et al. (2004) redefine the index to maximize 
the coverage of the field-of-view and retain the sensitiv- 
ity of the index towards changes in age, metallicity and 
abundance ratios. The new index name is Fe5270s. It 
measures the same spectral feature, but has a reduced 
spectral coverage in the red pseudo-continuum band. The 
new index can be converted to the original Lick/IDS sys- 
tem via the empirical form ula (Kuntschner et al. 2004, 
iFalcon-Barroso et alllioO^ : 

Fe5270 = 1.26 x Fe5270s -I- 0.12 

The 1(7 standard deviation of the above empirical calibra- 
tion is ±0.05 A for the Fe5270 index. More details on the 
derivation of the new index and its relation to the standard 
Lick/IDS index are given in Kuntschner et al. (2004). 

Having this in mind we measured Fe5270s and Mg6 
indices. The Mg6 index was measured using the Lick/IDS 
index definition, and all spectra were first broadened to 
the resolution of the Lick/IDS system. The Fe5270s index 
was later converted to index Fe5270 using above relation. 
Unfortunately, we were not able to determine the relevant 
offset to the Lick/IDS system, and correct for the sys- 
tematics, which come from differences in the continuum 
shape, because there are no reference stars in the HST 
archive observed by our and by the Lick/IDS instrumen- 
tal setup. The size of the corrections are probably similar 
to (or less than) our measurement errors. To first order, 
as well as for determining the relative trends in a galaxy. 
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Table 8. Line-strength indices measured in 0'.'55 x 0'.'2 aperture 



index 
slit 


NGC 4128 
cen neg pos 


NGC 4570 
cen neg pos 


NGC 4621 
cen neg pos 


NGC 5308 
cen —neg pos 


Mg& 
Fe5720 


5.0 ±0.2 5.0 ±0.4 

2.9 ±0.1 2.6 ±0.4 


4.5 ±0.2 2.9 ±0.3 3.7 ± 0.4 
3.1 ±0.1 1.7 ±0.3 1.5 ±0.4 


5.9 ±0.1 6.6 ±0.3 4.7 ±0.3 
3.7 ±0.1 3.4 ±0.3 3.4 ±0.3 


4.8 ±0.1 5.5 ±0.5 5.8 ±0.6 
3.1 ±0.1 4.1 ±0.4 3.5 ±0.5 


Mgb 
Fe5720 
B - V 


4.8 ±0.2 5.0 ±0.4 

2.8 ±0.1 2.6 ±0.4 

0.97 — 0.96 


4.4 ±0.2 2.9 ±0.3 3.5 ± 0.4 
3.1 ±0.1 1.7 ±0.3 1.4 ±0.4 


5.6 ±0.1 6.4 ±0.3 4.6 ±0.3 
3.6 ±0.1 3.4 ±0.3 3.3 ±0.3 
0.99 1.00 1.00 


4.7 ±0.1 5.3 ±0.5 5.7 ±0.5 
3.0 ±0.1 4.0 ±0.5 3.5 ±0.4 
0.96 0.94 0.96 



Notes - First two rows in the table present line-strengths (in A) corrected by the velocity dispersion measured using the whole 
spectral region. The second two rows present line-strengths (in A) corrected by the velocity dispersion measured excluding the 
Mgb region from the fit during extraction of kinematics. The last row presents B-V color measured at the actual positions of slits 
within the same slit-like aperture used for measuring line-strengths. The errors on the color values are estimated to be 0.05 mag. 



this is not very important, but has to be noted when com- 
paring with other studies. 

Broadening of the lines by the velocity dispersion 
weakens most of the lines and the index we mea- 
sure must be corrected for this effect. This can be 
achieved by determining an empirical correction factor 
C(cr)=index(0)/index(cr) for a star observed with the same 
instrumental setup. Index(O) is the index measured from 
the stellar spectrum, a is the velocity dispersion of the 
LOSVD with which the stellar spectrum is convolved and 
from which the index(cr) is measured. We used our un- 
broadened optimal template spectra to calculate the index 
at cr = and at the corresponding velocity dispersion, cr, 
of the galaxy spectrum. The proper correction factor C((t) 
was then applied to both measured indices. We used two 
approaches to extract kinematics and measure the veloc- 
ity dispersions (fitting the whole spectral region and ex- 
cluding Mgfe region from the fit). If the measured velocity 
dispersions differ, the velocity dispersion correction in the 
two cases will also be different. We noted the difference 
applying both corrections on the measured line-strengths. 

We measured the Mg6 and Fe5270 indices from each 
spectral bin used for kinematics. The corrected values and 
corresponding errors of the index are presented in the 
Tables A. 1-2 and shown in the Figs.|Hl-^2 The measured 
line-strengths for galaxies with higher S/N are relatively 
uniform with radius, rising towards the center, with dips 
in the case of the nuclei of NGC 4621 (Mgb) and NGC 5308 
(Fe5720). NGC 4128 does not show any trend, but rather a 
scatter of values, presumably due to the low S/N, while in 
case of NGC 5308 Mg5 line-strengths are slightly higher 
on one side of the galaxy. Generally, galaxies have high 
values of Mgb and Fe5270 indices. Detailed descriptions of 
spatially resolved line-strengths for all galaxies are given 
in Section^ 

The slits that were offset to the sides of the central 
slit do not have the required S/N to extract kinematics, 
but their summed spectra can be used to determine the 
indices on the positions of the galaxy outside the stellar 
disk. The binning of spectra is only useful up to the point 



at which summing more spatial elements does not simply 
add noise. We decided to use an aperture of 0'.'55 x 0'.'2 
(summing up 5 rows on each side of the central row in 
the spectral direction). With this approach the final side 
spectra used for the measurements of line-strengths had 
at least S/N « 10. In the case of NGC 4128, however, 
from one side spectra we were not able to extract any 
trustworthy measurement. Table |H| presents line- strengths 
corrected for velocity dispersions measured by fitting to 
the whole spectral range and excluding the Mg6 lines from 
the fit. As it can be seen, the differences between the lines 
are negligible and generally fall within the la error bars. 
Only in the case of the central slit of NGC 4621 which has 
the highest velocity dispersion as well as S/N, there is a 
significant difference. Adopting a conservative approach, 
we compared line-strengths from the second two rows in 
Table |S1 with the stellar population synthesis models in 
Fig.H 

We wish to compare the line-strengths measured on 
the disk with the line-strengths measured on the bulge 
using the three slit positions. The line-strength measure- 
ments on the summed spectra show similarly high values 
as the spatially resolved measurements, although with rel- 
atively lower values due to the smaller spatial resolution. 
The Mg5 index values in NGC 4621 are particularly high. 
Comparing with the literature w e find similar values for 
Mg6 and Fe 5270 index. Table 7 of'Traee r et~al\ l)l998|) list 
values of the same indices for NGC 4621 and NGC 4570 
(Mg& 5.50 and 4.65A, Fe5270 3.59 and 3.49 A), which, 
keeping in mind the unknown offset to the Lick s ystem 
and the lower spatial resolution of iTrager et all lll998^ 
data (aperture of 1" x 4"), are in good agreement with 
our findings. 

Age and metallicity have similar effects on the in- 
tegrated spectral energy distributions that we measure 
from unresolved sources due to a finely tuned conspiracy 
between age and metallicity variations (Worthev. ,199^ . 
Broad-band colors and many line-strength indices are de- 
generate with respect to age and metallicity. This makes 
the determination of the age and metallicities very diffi- 
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Fe5270 [A] 

Fig. 7. Age/met allicity diagnostic diagram (B-V color vs. 
Fe5270 index). Horizontal thick solid lines are lines of con- 
stant age [Gyr] and v ertical thin lines^ are lines of constant 
metallicity [Fe/H]) of lVazdekij l)l999() models. The size of 
symbols is related to the position of the slit: the smallest 
symbols are for the cen slits, intermediate for the pos slits, 
and the biggest for the neg slits. 



cult and ideally one would like to use two indices which 
can break this degeneracy. Usually, one or more Balmer 
lines {Hp, Hj, HS) are used as age indicators, and Mg6 
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dex values of our measured line-strengths also suggests 
the presence of non-solar abundances of elements. If not 
properly treated, over-abund ant indices can give wr ong 
age and metallicity estimates l|Kuntschner et al ] l200lh . A 
way around this issue is to define metal licity indicators 
which are insensitive to abundance ratios l|Gonzalejll993t 
iThomas et alll2003|) . The preferred indicator includes a 
combination of Mg&, Fe5335 and Fe5270 indices, where 
Fe5270 is the least sensitive to ch anges of [a/Fe] abun- 
dance ratios l|Thomas et aLlbOO^ . We were not able to 
construct such a metallicity indicator with the indices 
from our spectral range, and we chose to use the least 
sensitive Fe5270 index alone as a metallicity indicator. 
Since in our spectral range there are no age indicators, 
we decided to use a combination of broad-band B-V colors 
and Fe5270 index to construct an age/metallicity diagnos- 
tic diagram (Fig . [3. The model s presented by sohd lines 
are based on the lVazdeki3()l999(l single stellar population 
models: color values were obtained from A. Vazdekis web 
site'^, while we measured the Fe5270 index from the library 



\http:/ /www, iac. es/galeria/vazdekisj\ 



spectra broadening them to the Lick/IDS resolution. The 
combination of red colors and high metallicities puts the 
measured points on Fig. 13 on the top right of the model 
grid, indicating old stellar populations and a large spread 
of metallicities between the galaxies. 

4. Discussion 

In the two previous sections we presented the observa- 
tional results of HST program 8667. They include pho- 
tometric and spectroscopic observations of four galaxies 
with nuclear stellar disks. Here we analyze and discuss 
the observations. 



4.1. NGC 4128 

The most distant galaxy in the sample is NGC 4128 (36 
Mpc). It is an SO galaxy and it has not been detected in 
radio nor in IR. The isophotal parameters show that it is 
disky between 35 and 530 pc. On 1 kpc scale it has boxy 
parameters and on larger scales it becomes disky again. 

The color image shows a red nucleus. The values for 
Fe5270 index measured with an aperture bigger than the 
red nucleus are the smallest in the sample of galaxies. This 
combination puts the points on Fig. above the model 
grid. The difference in colors and line-strengths between 
the two slit positions are small and within errors indicate 
old stellar populations of ^^14 Gyr and metallicities be- 
tween [Fe/H] =-0.38 and solar. 

It is probable that the difference in the color between 
the nucleus and the rest of the galaxy, as well as the higher 
metallicity detected in the nucleus is connected with the 
unusual spatially resolved kinematic profiles (Fig.|Sl). The 
velocity dispersion is flat in the center. The velocity curve 
also shows an unusual flattening in the central 0'.'2, mea- 
surements being positive on the both sides of the galaxy 
nucleus. Outside this radii the galaxy rotates fast, as ex- 
pected for a disk galaxy. Kinematics extracted fltting to 
two different spectral regions are in a good agreement, 
confirming the results. Having in mind the boxiness in the 
central tens of arcseconds, the extracted kinematic indi- 
cates the existence of a small (~35 pc in diameter) core, 
kinematically distinct from the nuclear stellar disk. 

The B-V color profile on the last panel of Fig.Oshows 
a slightly shifted (~17 pc from the center) peak of the 
red nucleus. This supports the presence of a distinct com- 
ponent in the nucleus. On the other hand, the spatially 
resolved line-strengths do not follow this trend. The spec- 
tral observations of this galaxy have the smallest S/N ra- 
tio, and the significance of this discovery is just above 
la. Deeper exposures of high spatial resolution, prefer- 
ably with an integral-field unit to cover the 2D structure, 
are needed to confirm this result. 

4.2. NGC 4570 

NGC 4570 is a well-studied galaxy with HST. The main 
result from previous studies is that the inner region of 
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Fig. 8. Kinematic and line-strength profiles for 
NGC 4128. From top to bottom: mean velocity, ve- 
locity dispersion, Mgb index, Fe5270 index and B-V color 
profile. The color profile was extracted along the slit 
position averaging 0'.'2 perpendicular to the slit. The 
closed symbols represent measurement obtained by fitting 
the whole spectral region. The open symbols represent 
measurement by excluding the Mg& line from the fit. 



the galaxy was shaped under the influence of a weak bar 
I van den Bosch fc EmsellemlllQQsi) . The colors reveal no 
difference between the disk and the bulge, and a compar- 
ison with the stellar-population models indicate that the 
stars in the galaxy are of intermediate age, but the FOS 
spectral data gave a very high Hf] line-strength suggest- 
ing recent star formation (BJM98). One of the questions 
raised by these studies is whether all double-disk struc- 
tures are the result of bar-driven secular evolution. 

To the previous photometric and spectroscopic ob- 
servations we add new spatially resolved spectroscopic 
measurements with STIS (Fig. O. The probed region 
corresponds to the nuclear disk and inner 2". The ve- 
locity curve shows regular rotation peaking at ^ 0'.'15 
from the nucleus. The velocity dispersion steeply rises and 
peaks in the center. The kinematic profiles are similar to 
BJM98 ground based data, except the STIS data have a 
steeper velocity curve and higher velocity dispersion. In 
contrast, the FOS velocity dispersion from the same au- 
thors is about 50 km s^^ higher than STIS measurements; 
however, considering the error bars of both measurement 
(their error on sigma is « 30 km s~^) and the fact that our 
slit was significantly (for the width of the slit) offset from 
the galaxy nucleus, these measurements can be considered 
consistent with each other. 



Fig. 9. Kinematic and line-strength profiles for 
NGC 4570. From top to bottom: mean velocity, ve- 
locity dispersion, Mg6 index, Fe5270 index. The closed 
symbols represent measurement obtained by fitting 
the whole spectral region. The open symbols represent 
measurement by excluding the Mg5 line from the fit. 



In the case of NGC 4570, the central velocity dispersion 
is somewhat dependent on the spectral region used in the 
fit. Excluding Mg5 region systematically lowers the values 
by just over Icr, but increases the difference between this 
and BJM98 results. 

We also measure, within the errors, similar line- 
strengths to BMJ98, but with higher spatial resolution 
and we can give an estimate of the spatial changes in the 
indices. As can be seen from the Fig. both measured 
indices show flattening in the central 0'.'5 (~60 pc). At 
larger radii the metallicity drops. Also, the slits positioned 
on both side of the center measure the smallest metallicity 
(TablelHl) in the sample and the largest drop in values with 
respect to the center. This measurement shows that the 
nuclear disk consists of different stellar populations than 
the rest of the bulge, which is consistent with bar-driven 
evolution. 

4.3. NGC 4621 

The closest galaxy of the four is NGC 4621 (7 Mpc). It is 
also the only galaxy classified as an elliptical and is the 
only galaxy from the sample detected with IRAS (in the 
12 /im band). In the investigated range the galaxy is disky. 
The 64 coefficient steadily rises from the center to the dis- 
tance of 140 pc when it drops, but never reaching negative 
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Fig. 10. Kinematic and line-strength profiles for 
NGC 4621. From top to bottom: mean velocity, ve- 
locity dispersion, Mgb index, Fe5270 index and B-V color 
profile. The color profile was extracted along the slit 
position averaging 0'.'2 perpendicular to the slit. The 
closed symbols represent measurement obtained by fitting 
the whole spectral region. The open symbols represent 
measurement by excluding the Mg& line from the fit. 



values. With increasing radius it rises again, implying an 
outer disk. 

The color images reveal the most interesting features in 
the nucleus. The few central pixels are clearly much bluer 
than the rest of the bulge (TableEj). Another striking char- 
acteristic of the galaxy, mentioned in section 2.3.2, is the 
extended blue and red features visible on B-I (Fig.EJ and 
V-I images. The blue feature makes an angle of 15° with 
the major axis. Although the red feature spreads generally 
in the east-west direction (angle with major axis is 103°), 
it is not as clearly defined as the blue feature. We conclude, 
examining all color images, that the shape and the extent 
of the two features can be interpreted as a blue disk-like 
structure imbeded in the red bulge. The position angle of 
the blue feature is unexpected for an axisymmetric galaxy 
with a nuclear stellar disk. This significant structure per- 
haps can be explained considering the kinematics of this 
galaxy. 

The velocity and velocity dispersion panels in Fig. IIUI 
clearly show the existence of a kinematically decoupled 
core (KDC). This core was already detected by WEC02 
who observed the galaxy with integral field spectrograph 
OASIS mounted on CFHT and assisted by the PUEO 
adaptive optics system. Complementing their OASIS ob- 
servations, WEC02 also extracted kinematics from an 



archival STIS observations in the Ca triplet region, show- 
ing a peak in velocity dispersion C/.'OS from the center as 
well as confirming the KDC with the same spatial resolu- 
tion as in the data presented here. The difference between 
the WEC02 and our kinematic profiles is largest in the 
velocity dispersion profile, where our values lie systemat- 
ically above the WEC02 measurements. This difference is 
lower, but still present if we compare the kinematics ex- 
tracted excluding the Mgb line from the fit (~60 kms~^ 
for the velocity dispersion peak, but with large error-bars 
in both cases). This discrepancy could arise from the dif- 
ferent slit positions in the two studies, our slit being offset 
from the center and not covering the KDC uniformly. The 
counter-rotation of the KDC could lower the overall mea- 
sured dispersion if the slit is placed over its center, and, 
alternatively, if the slit misses the center of the KDC the 
measured velocity dispersion will be higher. 

The KDC on the WEC02 OASIS data is not ahgned 
with the major axis and it has a similar position angle 
as the blue feature on the color images presented here. 
Although the extent of the KDC is smaller (total of ^ 2" 
or 60 pc) than the blue feature on the B-I image, the 
existence of two structures could be connected as a result 
of the same formation process. 

Both Mg6 and Fe5270 line-strengths indices in 
NGC 4621 are the highest in the sample, also suggesting 
the over-abundance ratios of elemen ts simil ar to trends 
for g iant ellipticals ijKuntschned Il998 : Kuntschner et all 
I2OOII) . Our metallicity indicator, Fe5270, also has high 
values, with the central slit being slightly more metal- rich 
than the side slits as well as super-solar. Colors at the slit 
positions are red and the comparison with the stellar pop- 
ulation models indicates the age of the stars is between 10 
and 14 Gyr and the metallicity between solar and -1-0.2. 

Our spatially resolved measurements of the indices, 
show n in Fig. [ 101 are highe r than in previous studies 
(e.g. iKuntschner et alJ (|200ll) have central Mgb ^ 5.21 
with aperture of 3'.'4), but our aperture is much smaller 
(■~ 0'.'05) than that of any previous study, and the values 
outside the central arc-second approach the observed val- 
ues from the literature. The Mg6 index follows to some 
extent the changes in colors, showing a small dip in the 
center, while this can not be said for Fe5270 measure- 
ments. 

The existence of the KDC and the blue features in the 
red bulge of the NGC 4621 indicate two possible evolution- 
ary scenarios. The visible structures c ould be the result 
of a h ierarchical formation scheme (e.g. iKauffmann et all 
I1994I) involving a merger followed by a starburst where 
the KDC is the remnant of the ejected stars that later 
fell back in. These structures are r elatively long lived, 
having a relaxation time of ~1 Gyr (jBinnev fc Tremaind 
Il987l) : however, this is not long enough to explain the 
detected old age of the stars. Alternatively, the struc- 
ture could be produced by weak bar-driven evolution, 
as in the case of NGC 4570 l|van den Bosch fc EmsellemI 
'199^, where the observed double disk structure is the 
consequence of resonant frequencies in the galaxy, while 
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the blue feature and the KDC are the result of gas cap- 
tured on retrograde ('anomalous') or bits which are tilted 
with respect to the equatorial plane llPfenniger & Fricdli 
ll99lUFriedh fc Edrvlll99i lEmsellem & Arsenault. .19971 
Of course, a combination of both processes can also lead 
to the present situation. 

Distinguishing between the two scenarios is also dif- 
ficult because the galaxy has no obvious merger com- 
panion and it is nearly edge-on, making the detection 
of a weak bar more difficult. There are other cases of 
barred galaxies with similar properties to NGC 4621: i.e. 
an edge-on system with double-disk structure, unusual 
photometric and kinematic features and the difference in 
metallicity between the bulge and the disk. An example 
of a similar, although boxier, galaxy with a strong bar 
is NGC 7332. This galaxy was recently studied in de- 
tail by t^alcon-Barroso ct al. ( 2004). NGC 7332 is clas- 
sified as an SO galaxy and has a double disk structure 
l|Seifert fc Scorzalll99fil) . Examining their SAURON spectro- 
scopic observations, Falcon-Barroso et al. find a counter 
rotating stellar component within the central 250 pc. The 
galaxy also has complex gas morphology and the line- 
strength maps show it is young everywhere. The authors 
conclude that NGC 7332 is an SO galaxy with a bar viewed 
close to edge on. NGC 4621 and NGC 7332 are similar in 
their morphologies, and, although different in the stellar 
content, it is possible that NGC 4621 went through a sim- 
ilar formation process as NGC 7332. 

4.4. NGC 5308 

In many aspects NGC 5308 is different from the other 
galaxies in this study. Our photometry reveals the largest 
nuclear stellar disk in the sample of galaxies in this study. 
Unlike in the other galaxies, the nuclear disk of NGC 5308 
is very thin and bright. The diskiness parameter, b4, rises 
from the center and peaks at about 150 pc, dropping to 
zero at ~1 kpc and suggesting a distinction between the 
two disks. At large radii the galaxy again becomes disky. 

The stars in the disk of NGC 5308 rotate fast, reaching 
~100 km s~^ within 15 pc from the nucleus fFig. lll|l . The 
velocity dispersion has a peak of about 300 km s~^ in the 
center and is relatively flat in the inner 15 pc. This trend 
is also visible in the kinematics measured excluding the 
Mg6 line from the fit, although the right hand side of the 
plot shows considerably lower velocity dispersion values. 
This is reflected in the panel with Mgb values, which are 
slightly higher on the right hand of the plot. The Mg6 
and Fe5270 index values have opposite trends. The small 
variations in the B-V color along the slit, including the 
sudden blue dip in the center, are followed by the line- 
strength measurements. 

There is a big difference between the line-strengths 
measured on the different slit positions (Table [S] and 
Fig. O. The nucleus, being just below solar metallicity, 
is more metal-poor than the bulge which has a non-solar 
abundance ratio of elements and the highest metallicity 



in the sample. The nucleus and the investigated part of 
the bulge also have different colors, with the center being 
redder. Comparing these results with the stellar popula- 
tion synthesis models reveals an old stellar population in 
the nucleus (14 Gyr), while the colors of the bulge suggest 
intermediate age stellar component (~ 5-10 Gyr). 

The color difference between the bulge and the nucleus 
is measured because the side slits were positioned in the 
region of the minor axis blue feature, especially visible 
on B-I image. This region is also special for its high metal 
content. There are no hints of specific morphological struc- 
tures (such as a polar ring) along the minor axis and the 
question is whether the rest of the bulge, especially the 
red feature, and the nuclear stellar disk share the same 
metallicity. Unfortunately our spectral observations did 
not cover the necessary areas and we can only speculate 
on the processes that created the observed structures. 

If the metallicity of the minor axis blue feature and 
the nuclear disk are the same, the formation of these two 
structures has to be connected to the same formation sce- 
nario, most probably involving a transp ortation of gas to 
the ce nter, perhaps by a bar, for which ISeifert fc Scorzal 
||1996^ found evidence in NGC 5308. If the metallicity of 
the nuclear disk is lower than the metallicity of the minor 
axis blue feature, but similar to the measured metallicity 
of the nucleus, the two blue features in NGC 5308 were 
not created from the same infalling material, but they still 
can be from the same epoch. If the metallicity of the mi- 
nor axis blue feature is equal to the rest of the bulge, 
which is redder and therefore older, then the younger stel- 
lar population in the blue feature must have been induced 
by an internal process, perhaps ionization from the ra- 
diation generated by the central black hole, which was 
turned on with the infall of the material that made the 
blue disk , and late r turned off with the stabilization of 
the disk ijLoeb fc R asio 1994). 

In section 2.3.2 we showed that the major axis blue fea- 
ture in NGC 5308 corresponds to the nuclear stellar disk. 
It is about 30 pc thick which strongly suggest the galaxy is 
viewed very close to edge-on. Comparing wit h the vertical 
scalelength of 34 edge-on spirals presented bv'Kre gel et alJ 
I2OO2), which are between 0.2 and 1.4 kpc thick, the disk 
in NGC 5308 is a remarkably thin disk. Note that our 
estimate of the disk thickness can only be approximately 
compared with the vertical scalelength measurements of 
Kregel et al. It is also not possible to say much about 
the thickness of the other nuclear disks, due to their in- 
clination (not as edge-on as NGC 5308). Whether this 
nuclear disk resembles the disks from the group of "super- 
thin" galaxies, like UG C 7321 or IC 5249 l lVIatthe ws et all 
»1999£ .Matthews 2000: Ivan der Kr^iit et al. 2001*), is an 
open question. The sizes of the nuclear disk in NGC 5308 
and known "super-thin" disks are quite different as well 
as the surrounding environment (stellar bulges and dark 
matter halos respectively). A proper way to compare the 
disks is to measure the radial and vertical sizes in a con- 
sistent manner, which is beyond the scope of this paper. 
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Fig. 11. Kinematic and line-strength profiles for 
NGC 5308. From top to bottom: mean velocity, ve- 
locity dispersion, Mgb index, Fe5270 index and B-V color 
profile. The color profile was extracted along the slit 
position averaging 0'.'2 perpendicular to the slit. The 
closed symbols represent measurement obtained by fitting 
the whole spectral region. The open symbols represent 
measurement by excluding the Mg& line from the fit. 

The color and metallicity of the disk in NGC 5308 
suggest the disk could be made of a younger and more 
metal-poor stellar population than the rest of the galaxy, 
implying it formed at a different epoch from accreted ma- 
terial. 

5. Conclusions 

We have presented photometric and spectroscopic ob- 
servations of four nearby early-type galaxies with nu- 
clear stellar disks (NGC 4128, NGC 4570, NGC 4621, 
NGC 5308). The observations consist of high resolution 
images with WFPC2 using the F450W and F555W fil- 
ters, and STIS high resolution spectra through the 52x0.2 
long-slit with the G430M prism. 

The photometric analysis reveals similarities and dif- 
ferences between the galaxies. Nuclear stellar disks are 
clearly visible and are photometrically disconnected from 
the large scale disks. NGC 4128 shows boxy isophotes on 
the inner and outer edge of the nuclear disk. NGC 4621, 
the only E galaxy in the sample, is everywhere disky, while 
NGC 5308 has a razor-thin 30 pc) disk. 

Color images reveal interesting and unexpected struc- 
tures. NGC 4128 has a red nucleus, while NGC 4621 has 
a blue nucleus. Prominent color features are visible on all 
galaxies. The blue feature in NGC 4128 is analyzed inlBl 



and is the signature of a transient event. NGC 4621 has a 
blue feature at an angle of 15° with the major axis on top 
of a red bulge. It is likely connected to the KDC discovered 
by WEC02. The nuclear stellar disk in NGC 5308 is as- 
sociated with the razor thin blue feature along the major 
axis. NGC 5308 has another blue feature along the minor 
axis. The colors of all three galaxies indicate old stellar 
populations except for the bulge of NGC 5308 where the 
combination of slightly less red colors and high metallicity 
lowers the age of the stellar populations. 

The high resolution spectroscopy was obtained at three 
positions on each galaxy. One slit was positioned on the 
nuclear stellar disk, with the PA equal to the major axis 
PA. Two additional slits were positioned on both sides of 
the central slit, ^ 0'.'3 away from the disk covering the 
bulge. The central slits were used to extract spatially re- 
solved kinematics and line-strengths. The S/N permitted 
the extraction of the mean stellar velocity and the velocity 
dispersion, as well as the measurement of line-strengths. 
The kinematics will be used in a separate paper to esti- 
mate the black hole masses in the centers of the galaxies. 

Considering the shape of the spatially resolved kine- 
matic curves, the four galaxies could be sorted in 
two groups: fast and kincmatically disturbed rotators. 
NGC 4570 and NGC 5308 belong to the first group. Their 
rotation curves show clear signature of the stellar disks. 
The rotation curves of NGC 4128 and NGC 4621 are much 
more complicated. In the case of NGC 4621 the unusual 
mean velocity and velocity dispersion curves are consistent 
with the known KDC (WEC02) in the nucleus. Although 
based on a la detection, we report the discovery of a sim- 
ilar kincmatically distinct core in the case of NGC 4128. 

Spatially resolved line-strength measurements along 
the disk indicate that all four galaxies are more metal-rich 
in the inner 0'.'5 than outside this radius. Both measured 
indices (Mg6 and Fe5270) increase towards the center, ex- 
cept in the case of NGC 5308 where Fe5270 has an op- 
posite trend to Mg6 index. Non-solar abundance ratios of 
[Mg/Fe], hinted by results of the extraction of kinemat- 
ics, are present in NGC 4570 and NGC 4621, and to some 
extent also in NGC 5308. 

The slits positioned on the bulges had low S/N and 
no spatially resolved kinematics were extracted. However, 
by binning the spectra, it was possible to measure the 
line-strengths at one position and compare them to the 
values of the nuclei. The objects show various structure: 
NGC 4128 has similar metallicities at the different slit po- 
sitions, NGC 4570 and NGC 4621 have higher metaUicity 
in the nucleus, while NGC 5308 in the bulge. Generally, 
the galaxies show a spread in metallicity from sub- to 
super-solar. 

This study shows the diversity within this class of ob- 
jects, but also emphasizes the similarities in the photom- 
etry and kinematics. The red color gradient in the nu- 
clei of NGC 4128 and the blue features in NGC 4621 
and NGC 5308 suggest the existence of different stel- 
lar populations on small scales ('^100-500 pc). The in- 
vestigated galaxies were chosen as galaxies with specific 
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nuclear morphologies: nuclear stellar disks. However, ex- 
cept in the case of NGC 5308 the colors of the disks are 
not much differe n t from the bulge, as previously noted 
by ICarollo et all lll997bl) . The existence of other color 
features is a surprise. In two galaxies (NGC 4128 and 
NGC 4621) these color features are followed by the ex- 
istence of a KDC. The other two galaxies do not show 
any peculiarities in their kinematics. Also, if the KDC in 
NGC 4621 is connected to the misaligned blue fe ature, 
we can conclude, similar to ICarollo et al.l l(l997al) . that 
KDCs are not kinematic counterparts of the nuclear stel- 
lar disks. This gives credit to the complexity of forma- 
tion scenarios that demands a separate study per galaxy, 
b ut there are a few most likely fr ameworks, outlined also 
in lScorza fc van den BoschI l)l998l) . in which the processes 
responsible for the observed structures operate. 

The formation of nuclear disks, rings and double 
disk structures in early-type galaxies can be explained 
through secular evolution dri ven by w eak bars as shown 
by Ern scllem ct, al. ( 1996) and ^van den B osch & EmscUcni. 
l|l998j) in the cases of M 104 and NGC 4570, respec- 
tively. This mechanism, through the evolution of the 
bar, explains the double-disk morphology. Support for 
this scenario comes from the fact tha t SO galaxies have 
high line-strengths ((Fisher et al.lll99'6l) and there are evi- 
dences of embedded ba rs in early-type gala xies (e.g. M104, 
NGC 4570, NGC 7332. IScorza et alJll998l) . This model of 
bar-driven evolution is consistent with the observations 
in the presented galaxies, even in the cases of the galaxies 
with KDCs, such as NGC 4621, but also NGC 4128, which 
is additionally boxy and presents an interesting case. The 
time varying triaxial potential of the bars offers exotic 
orbits that could explain the existence of kinematic and 
photometric features. In this scenario, the KDCs are cre- 
ated from enriched material transported inwards (perhaps 
even gas acquired through a merger) , which gets frozen on 
retrograde orbits tilted with the respect to the equatorial 
plane. 

Other possibilities involve a merger scenario (capture 
of gas that settles in the principle plane forming stars, 
and/or makes tidal inflow s that create KDCs), or growth 
of a central black hole l|Loeb fc Rasid Il994|) . A black 
hole stabilizes the disk and within this scenario a con- 
nection to quasars can be made by stopping the fuel- 
ing of the central engines with the formation of a sta- 
ble disk. None of the previously investigated nuclear stel- 
lar disk galaxies has an acti ve nucleus, although they 
do harbor IO^m^ black holes ( Kormendvet al.1 Il996blat 
ICretton fc van den Boschll999|) . This makes them descen- 
dants of quasars that spent their fuel (there is not much 
dust or gas in most of these galaxies), or quasars that, 
through dynamical evolution, turned off the central engine 
(stabilization of the disk due to the growth of the black 
hole, disappearance of bars that transport the material to 
the center). 

Nuclear disks are easier to find in edge-on systems; 
however, the influence of weak bars is correspondingly 



more difficult to ascertain. Detailed spectroscopic stud- 
ies with two dimensional coverage of the major features 
(nuclei, stellar disks, KDC, photometric features) are nec- 
essary to chose between the present formation scenarios. 
A careful investigation of the two-dimensional kinematic 
properties and their connection to the distribution of line- 
strengths (metal content and age of stellar populations) 
can offer decisive tools to deduce the nature and nurture 
of galaxies with nuclear stellar disks. 
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Appendix A: Extracted kinematics 



Table A.l. Measured kinematics and line-strengths for 
observed galaxies using full spectral region in the fit 



radius 


V 


5V 


(J 


5g 


Mg6 


5Mgb 


Fe5270 


(5Fe5270 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 4128 


-0.95 


2323. 


29. 


168. 


46. 


3.6 


0.9 


4.1 


0.5 


-0.50 


2359. 


19. 


258. 


27. 


5.3 


0.7 


4.3 


0.4 


-0.30 


2438. 


24. 


283. 


27. 


4.3 


0.6 


3.6 


0.4 


-0.15 


2370. 


20. 


206. 


24. 


5.1 


0.6 


2.5 


0.4 


-0.10 


2475. 


20. 


253. 


31. 


4.1 


0.6 


3.5 


0.3 


-0.05 


2489. 


12. 


229. 


15. 


4.6 


0.5 


2.9 


0.3 


0.00 


2474. 


13. 


230. 


14. 


5.6 


0.5 


2.8 


0.3 


0.05 


2498. 


14. 


216. 


18. 


5.6 


0.5 


3.4 


0.3 


0.15 


2568. 


14. 


191. 


11. 


4.6 


0.5 


3.3 


0.3 


0.25 


2618. 


12. 


150. 


13. 


4.7 


0.7 


1.7 


0.4 


0.50 


2649. 


18. 


203. 


27. 


5.7 


0.8 


1.5 


0.5 


1.25 


2691. 


20. 


91. 


30. 


4.1 


1.0 


2.2 


0.7 


NGC 4570 


-0.70 


1827. 


13. 


138. 


13. 


2.0 


0.5 


2.4 


0.4 


-0.45 


1803. 


16. 


106. 


18. 


2.0 


0.6 


2.2 


0.3 


-0.25 


1853. 


9. 


136. 


9. 


4.1 


0.4 


2.6 


0.3 


-0.15 


1856. 


9. 


155. 


14. 


4.3 


0.4 


3.3 


0.3 


-0.05 


1877. 


8. 


169. 


8. 


4.8 


0.4 


3.4 


0.3 


0.00 


1920. 


8. 


202. 


12. 


4.8 


0.4 


3.6 


0.3 


0.05 


1941. 


8. 


167. 


9. 


4.5 


0.4 


2.7 


0.3 


0.10 


1950. 


7. 


145. 


9. 


4.7 


0.4 


3.3 


0.3 


0.15 


1979. 


7. 


123. 


10. 


4.3 


0.4 


3.3 


0.3 


0.30 


1966. 


13. 


176. 


17. 


4.5 


0.5 


3.0 


0.3 


0.50 


1996. 


10. 


133. 
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3.3 


0.5 


2.2 


0.3 


0.80 


1956. 


16. 
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3.2 


0.6 


2.1 


0.4 


NGC 4621 
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335. 
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0.5 
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-0.50 


572. 
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306. 
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546. 
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315. 
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308. 
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3.0 
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3.4 


0.2 


-0.05 


631. 


13. 


339. 


13. 


5.4 


0.3 


4.0 


0.2 


0.00 


594. 


13. 


356. 


15. 


5.8 


0.3 


4.0 


0.2 


0.05 


603. 


17. 


420. 


29. 


6.4 


0.3 


3.8 


0.2 


0.10 


588. 


24. 


395. 


34. 


6.2 


0.4 


3.8 


0.2 


0.20 


629. 


13. 


341. 


21. 


5.7 


0.3 


4.0 


0.2 


0.35 


605. 


12. 


283. 


10. 
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0.4 


4.0 
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0.65 


552. 


14. 


283. 


14. 


5.0 


0.4 


3.2 


0.3 


NGC 5308 


-0.90 
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-0.50 


2241. 


20. 
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26. 
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0.6 
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0.3 


-0.30 
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11. 


221. 


15. 
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0.4 
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0.2 


-0.15 


2259. 


13. 
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12. 


4.1 
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3.0 


0.3 


-0.10 


2258. 


12. 


279. 


17. 
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3.4 


0.2 


-0.05 


2179. 


14. 
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21. 
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0.3 


2.6 


0.2 


0.00 
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16. 
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0.3 


2.9 


0.2 


0.05 


2007. 
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17. 
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0.3 


2.9 
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0.10 


2003. 


11. 
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23. 
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2.6 


0.3 


0.20 
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9. 


175. 


12. 
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0.4 


3.0 
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0.30 
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8. 


147. 


10. 
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0.5 


3.3 
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10. 
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11. 


4.4 
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3.3 


0.3 



Notes - Col. (1): radius of measurement (arcsec); Col. (2): 
mean velocity (kms~^); Col. (3) mean velocity error (kms~^); 
Col. (4): velocity dispersion (kms"'^); Col. (5): velocity disper- 
sion error (kms^^); Col. (6): Mg6 index (A); Col. (7) Mg& in- 
dex error (A); Col. (8) Fe5270 index (A); Col. (9) Fe5270 index 
error (A) 



Table A. 2. Measured kinematics and line-strengths for 
observed galaxies excluding Mg6 region from the fit 



radius 


V 


5V 




Sa 


Mgfo 5Mg6 Fe5270 5Fe5270 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 4128 


-0.95 


2294. 


31. 


258. 


30. 


3.9 


0.8 


4.4 


0.5 


-0.50 


2378. 


25. 


275. 


34. 


5.4 


0.6 


4.3 


0.4 
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2431. 


35. 
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30. 


228. 
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258. 
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229. 


24. 


4.6 
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2495. 


19. 


223. 


25. 


5.5 


0.5 
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0.3 


0.05 


2520. 


19. 


223. 


29. 


5.6 


0.5 


3.4 


0.3 


0.15 


2603. 


15. 


174. 


12. 


4.4 


0.5 


3.2 


0.3 


0.25 


2674. 


15. 


118. 


15. 


4.4 


0.7 


1.6 


0.4 


0.50 


2657. 


24. 


179. 


24. 


5.6 


0.7 


1.4 


0.6 


1.25 


2688. 


23. 


74. 


33. 


4.0 


1.0 


2.2 


0.7 


NGC 4570 


-0.70 


1805. 


14. 


138. 


14. 


2.0 


0.5 


2.4 


0.4 


-0.45 


1815. 


21. 


131. 


23. 


2.0 


0.6 


2.2 


0.4 


-0.25 


1853. 


8. 


133. 


10. 


4.1 


0.4 


2.6 


0.3 


-0.15 


1872. 


9. 


142. 


15. 


4.2 


0.4 


3.2 


0.2 


-0.05 


1889. 


9. 


149. 


10. 


4.8 


0.4 


3.4 


0.3 


0.00 


1921. 


11. 


176. 


18. 


4.7 


0.5 


3.6 


0.3 


0.05 


1942. 


9. 


146. 


9. 


4.4 


0.4 


2.7 


0.3 


0.10 


1957. 


8. 


118. 


8. 


4.7 


0.4 


3.3 


0.3 


0.15 


1993. 


8. 


93. 


6. 


4.2 


0.4 


3.3 


0.3 


0.30 


1985. 


14. 


142. 


16. 


4.4 


0.5 


3.0 


0.3 


0.50 


1995. 


11. 


123. 


10. 


3.3 


0.5 


2.2 


0.3 


0.80 


1988. 


17. 


111. 


25. 


3.1 


0.6 


2.2 


0.4 


NGC 4621 


-0.95 


531. 


27. 


306. 


29. 


4.3 


0.5 


3.6 


0.3 


-0.50 


569. 


19. 


203. 


20. 


4.9 


0.4 


3.2 


0.3 


-0.30 


562. 


21. 


280. 


19. 


5.6 


0.4 


3.0 


0.3 


-0.15 


612. 


18. 


278. 


13. 


5.6 


0.3 


2.9 


0.2 


-0.10 


622. 


16. 


316. 


15. 


6.0 


0.3 


3.3 


0.2 


-0.05 


663. 


15. 


323. 


14. 


5.1 


0.3 


4.0 


0.2 


0.00 


608. 


17. 


307. 


17. 


5.5 


0.3 


3.8 


0.2 


0.05 


589. 


26. 


379. 


30. 


6.1 


0.3 


3.7 


0.2 


0.10 


592. 


39. 


324. 


51. 


5.7 


0.4 


3.5 


0.2 


0.20 


659. 


17. 


288. 


20. 


5.3 


0.3 


3.8 


0.2 


0.35 


618. 


14. 


293. 


12. 


5.5 


0.4 


4.1 


0.2 


0.65 


531. 


24. 


286. 


20. 


5.0 


0.4 


3.2 


0.3 


NGC 5308 


-0.90 


2333. 


24. 


110. 


37. 


3.9 


0.7 


2.2 


0.4 


-0.50 


2188. 


22. 


266. 


21. 


3.7 


0.5 


4.1 


0.3 


-0.30 


2256. 


14. 


244. 


19. 


4.6 


0.4 


4.1 


0.2 


-0.15 


2270. 


12. 


181. 


10. 


4.0 


0.4 


2.9 


0.3 


-0.10 


2261. 


16. 


291. 


26. 


4.7 


0.3 


3.5 


0.2 


-0.05 


2173. 


22. 


290. 


32. 


4.9 


0.3 


2.5 


0.2 


0.00 


2122. 


21. 


305. 


20. 


4.8 


0.3 


2.9 


0.2 


0.05 


2001. 


14. 


221. 


42. 


5.2 


0.3 


2.8 


0.2 


0.10 


2042. 


13. 


96. 


25. 


4.8 


0.4 


2.5 


0.3 


0.20 


1954. 


9. 


123. 


13. 


4.9 


0.4 


2.9 


0.3 


0.30 


1965. 


9. 


136. 


11. 


4.8 


0.5 


3.3 


0.3 


0.55 


1955. 


10. 


115. 


16. 


4.3 


0.5 


3.2 


0.3 



Notes - same as in Table rO 



22 



Krajnovic & JafFc: HST observations of nuclear stellar disks 





Fig. B.l. Zoom-in on the WFPC2 images of the nu- 
cleus of the NGC4128, in F702W (left) and F555W band 
(right). Both images are oriented as on Fig. El (north is 
down and cast to the right). Horizontal line is 1" long. 
The left image was taken almost two years before the right 
one. The offset of the peak next to the center on the right 
image is 0'.'143 west and 0'.'318 north from the center of 
the galaxy. 

Appendix B: A transient in NGC 4128 

Examining the HST - WFPC2 images of the nearby SO 
galaxy NGC 4128, taken on May 14 2001 in the F450W 
and F555W filters, a peculiar source of light was discov- 
ered very close to the center of the galaxy fFig. lB.lT) . The 
source is offset by 0'.'14 west and 0'.'32 north of the galaxy 
center on the WFPC2 images. In the HST archive there 
are images of the galaxy taken on March 17 1999, with fil- 
ter F702W, observed under FID 6357 (PI Jaffc), which do 
not show any additional light source next to the nucleus 
of the galaxy. This clearly indicates the appearance of a 
transient object in our data. In this Appendix we try to de- 
duce the origin of the transient. Section rB.ll presents a new 
method of extraction of the light contribution from the 
background galaxy and the photometric analysis, followed 
by discussion on the nature of the object in section Ir!21 

B.l. Extraction of light and results 

The general properties of NGC 4128 were described in the 
main text. Here we add that the du st in the galaxy has 
not been detected ijTran et al.ll200l|) . This ensures there 
is no significant extinction although the object is almost 
in the center of the galaxy. The galactic component of the 
reddening in the direction of the NGC 4128 is es timated 
to be E{B -V)^ 0.02 mag llSchlegel et al.lll998l) . 

To determine the nature of the transient, it is neces- 
sary to accurately do photometry on the observed light. 
The first required step is to remove the contribution of the 
galaxy light. Since there are no pre-transient observations 
of the galaxy in the two filters, it was not possible to use 
the standard technique of subtracting a suitable image of 
the galaxy taken with th e same in strumental set up (e.g. 
iFilippenko et al.lil986: Hamuv et al., 1994) . To overcome 
this problem a Multi-Gaussian Ex pansion (MGE) mode l 
of the galaxy was constructed (e.g.'Emsellem et al."'l994'). 
The method and software provided bv .Cappellari (,2002) 




Fig. B.2. Contour maps of WFPC2 F450W (left) and 
F555W (right) images of NGC4128. The images show the 
inner 7'.'5. Superposed on the two images are the con- 
tours of the MGE surface brightness, convolved with the 
WFPC2 PSFs for the corresponding filters. The transient 
was excluded from the fit, but it is shown in the images 
for illustration. 



were used. The MGE models of the galaxies were created 
using only the PCI chip and masking the region of 
the transient by a circular mask with radius of 5 pixel 
(Fig. |R2|| . For both MGE models, the total number 
of used Gaussians was increased until the minimum 

stopped decreasing. This approach yielded the total 
number of 9 Gaussians and RMS error of about 2%. 

The obtained models were subtracted form the original 
images. Figure IbTsI shows a horizontal cut through the ob- 
served image at the position of the transient and the corre- 
sponding MGE model. The light of the transient is clearly 
distinguishable from the rest of the galaxy. Subtracting 
the MGE model of the galaxy from the observations it is 
possible to recover the transient light contribution. The 
residuals are symmetric and can be used to calculate the 
magnitude of the transient. Photometric measurements 
were conducted using the IRAF task PHO T. Magnitud e 
zero points for filters were obtained from iDolohinl l|200o|) . 
The values calculated in an aperture of 0"5 were corrected 
to infinite aperture by adding 0.1 mag. The correction for 
geometrical distortion and the CTE correction were not 
performed as their infiuence are of second order, com- 
pared to the uncertainty introduced by subtracting the 
model galaxy from the original image. The final transient 
parameters, namely its position, apparent magnitude and 
absolute magnitude in Johnson-Cousins system, are given 
in the Table EH 

B.2. Discussion and conclusions 

Lacking any spectra of the detected object as well as a 
longer period of observation, it is difficult to determine the 
object's true nature. Since it is not visible on the F702W 
band image while it is quite bright on the F450W and 
F555W band images (taken 22 months later), it is most 
natural to conclude the object is a transient. There are a 
few possibilities such as a solar system object, a nova or 
a supernova. 
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Table B.l. Summary of transient measurements. 
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filter 


ra 


dec 


S X 




counts 


m 


0"m 


M 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F450W 


f2:08:34.84 


68:46:57.15 


0.142 


0.315 


20940.09 


19.06 


0.05 


-13.7 


F555W 


f2:08:34.84 


68:46:57.25 


0.143 


0.318 


18099.38 


19.17 


0.05 


-13.6 



(1): name of the filter; (2) and (3): position of the transient measured in (h,m,s) and (deg,arcmin,arcsec) on WFPC2 (J2000); 
(4) and (5): the offset in arcsec of the transient from the intensity peak of the galaxy due west and north respectively; (6) 
total counts inside aperture of 0'.'5; (7) magnitude of the transient in Johnson-Cousins system; (8) the estimated error on the 
apparent magnitude; (9) absolute magnitude using the distance to NGC 4128 from Table [T] 





Fig. B.3. Horizontal cuts through the galaxy on the po- 
sition of the SN. The thick line is the MGE model and 
the thin line is the observed light. The lower plots show 
residuals obtained subtracting the models from the obser- 
vations. The transient is clearly separable from the rest of 
the galaxy by a MGE model. 

The duration of all four independent WFPC2 expo- 
sures is about 40 minutes. A solar system object would 
show a noticeable movement between the different expo- 
sures, except if moving directly toward an observer. If the 
object is at the distance of the Kuiper belt it would move 
approximately a few arc seconds during the observations. 
However, the relative position of the object (the distance 
between the object and the center of the galaxy) changes 
by less than one pixel (< 0'.'0455), the movement being less 
than measurement error. It seems reasonable to conclude 
the object is of extra-Solar system origin. 

Novae are known to have a maximum absolute 
magnitude less bright then -9 mag iCoheaiaS^. With 
the measured absolute magnitude of -13.6 (V band), 
this is ruled out. Similarly, because it is very faint, 
the observed transient is also not likely a nova in our 
galaxy (a compilation of novae light curves is given in 
Ivan den Bergh fc Younge^ ijlQSTfl ). The observed magni- 
tude leaves the possibility that the object in NGC 4128 is 



a supernova. The type of supernova is defined according 
to its spectrum, and here again it is hard to establish any- 
thing specific. However, the host galaxy is an SO galaxy, 
and since supernovae of t ype II do not occur in early- 
type galaxies (Table 3. in ICappellaro et allllQOTbl) . the 
transient is probably a supernova of Type la. Comparing 
the measured absolute magnitu de to absolute magnitude s 
of supernova 1994D or 1992 A jCapDellaro et alJll997ah . 
which were also SO galaxies, suggests that the supernova 
was observed about 200 days after the explosion. The 
B — V color of the supernova is -0.11, which is bluer 
then the expected color for a type la SN at B maximum. 
On the other had, it is consistent with the B — V light 
curve at the later time (around 200 days) of the above 
mentioned reference supernovae (ISalvo et alJl200ll) . 

Private communication with the supernova-survey 
groups (D. Green from lAU Circulars and W. Li from 
LOTOSS survey) did not confirm the existence of the 
supernova. They checked images from January to June 
2001, and after galaxy subtraction there was no signature 
of any transient object. However, the proximity to the 
nucleus of NGC 4128 is a possible reason why ground- 
based automated search program might have overlooked 
the supernova. 

If the source of the light is a Type la supernova, 
it will be the first supernova discovered in NGC 4128 
ijBarbon et alJll999l) "^. Unfortunately, without other high 
resolution available data it is impossible to determine the 
true nature of the object. Also, the estimate of the timing 
of the supernova in NGC 4128 is highly uncertain, since 
it depends on precise classification and calibration. Type 
la supernovae are known to have different absolu te mag- 
nitude light curves (e.g. ICappellaro et al.lll997a|^ . If the 
transient in the NGC 4128 is not a supernova of type la 
200 days after explosion then its nature remains unknown 
(a flaring second black hole, an asteroid on the trajectory 
to Earth). 



http : //web .pd. astro . it/ supern/ 



